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Abstract 

This note shows that for water droplets at weakly absorbing wavelengths, the ratio 
wox(r)/a)oA(ro) of two single scattering albedo spectra, (jooa(V) and cooaOV), is a linear 
function of o)oa(/). The slope and intercept of the linear function are wavelength 
independent and sum to unity. This relationship allows for a representation of any single 
scattering albedo cooA(r) via one known spectrum gooa(/o). The note provides a simple 
physical explanation of the discovered relationship. In addition to water droplets, similar 
linear relationships were found for the single scattering albedo of non-spherical ice 
crystals. 
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1. Introduction and the main result 

The single scattering albedo in atmospheric radiative transfer is the ratio of the scattering 
coefficient to the total extinction coefficient (e.g., Bohren and Clothiaux, 2006, p. 257). It is 
equal to unity if all extinction is due to scattering (e.g., Rayleigh molecular scattering); 
conversely, it is equal to zero if all extinction is due to absorption (e.g., gaseous absorption) 
For spherical particles with a given refractive index single scattering albedo is calculated 
using Mie theory. For cloud water droplets both the scattering and absorption coefficients, 
thus the single scattering albedo, are functions of wavelength and droplet size. If the 
distribution of droplet sizes is given, both the coefficients and the single scattering albedo 
can be determined as functions of effective radius provided that the width of droplet size 
distribution is constant. 

In general, the single scattering albedo ooo, as a function of effective radius r and 
wavelength A, is defined as, 
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where ct s a and cr a A are the scattering and absorption coefficients, respectively. To compute 
the single scattering albedo a Mie code (Bohren and Huffman, 1983) was run for four 



wavelengths (0.86, 1.65, 2.13 and 3.75 pm) and for r from 5 to 50 gm (a gamma size 
distribution with effective variance 0.1 gm (Hansen and Travis, 1974) was assumed). 
When a ratio a)oA(r)/a>oA(ro) with ro=5 pm was plotted against oooA(r) as a function of the 
four wavelengths A, a remarkable linear relationship was discovered: 
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Here slope pi and intercept p 2 are wavelength independent parameters. Since Q)oA(r)=l at 
A=0.86 gm for all effective radii r, the sum of pi and p 2 is equal to 1, i.e. (2) can be rewritten 
as 
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Figure 1 illustrates Eq. (3) for r - 10, 20, 30 and 40 gm and ro=5 gm. In the next section we 
provide a simple physical explanation of the spectral invariant relationship (3) and justify 
the results illustrated in Fig. 1. 


2. Justification of the main results for water droplets 

We start with the definition of the single scattering albedo (1). For the visible and near- 
infrared spectral regions, the scattering coefficient does notvary much with either 
wavelength or droplet size; we can assume it to be constant, i.e., 
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Thus Eq. (1) can be rewritten as 
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and the ratio of co OJL for two effective radii r and ro (ro < r) will be 

“oWo) + „x(r) . (6) 

Following Twomey and Bohren (1980), for weakly absorbing water droplets 

l-o)oA(r) ^ kxr (7) 

where k\ is the bulk absorption coefficient (4tt multiplied by the ratio of the imaginary part 
of the refractive index to wavelength). For wavelengths A between 0.2 and 2.5 gm, the 
absorption coefficient 10‘ 4 </a<10 2 cm' 1 (Bohren and Clothiaux, Fig. 2.25, pg. 113). 



Next, since 
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it follows from Eq. (7) that 
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and, neglecting the second order terms o[(far] 2 ], we get 

r 


a al {r)~0,-+ r = 0 l [k l r + {kj) +(k t r) +...] = ff/,r 


Let us now define a wavelength independent parameter p as 

p = 1-ro/r. 


It follows that 


r 0 = r(l-p). 


We now rewrite the ratio (6) as 
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Applying approximation (10] and relationship (12] to the second term in (13] gives 
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Equation (14] completes the "proof' of (3] with p defined by (11]. Equation (3] is an 
approximation but it can be good enough for many purposes. Rearranging Eq. (3] leads to 
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Equation [15] states that the single scattering albedo co 0X (r] can be expressed via one 
known spectrum c%(ro]. Indeed, if one knows r and ro, the slope p=p[r,ro) can be 
approximated by Eq. (1 1). However if r and ro are unknown, p can be obtained directly 
from Eq. (3). For example, if one knows the two single scattering albedos gooa(/] and 
u)ox(ro) at least at one weakly absorbing wavelength, say, A=2.1 pm, then p can be 
determined from Eq. (3) as the slope between two points at A=2.1 pm and at any visible 
wavelength since moA(r]=(jOoA(ro]=l there. Knowing p, Eq. (15] provides c%(r] for any 
wavelength A if c%(ro] is available. 

How well does Eq. [11] approximate the slope p as a function of r and ro? Figure 2 
illustrates two slopes: p from the Mie calculations and its approximation p a p P r=l-r/ro as a 
function of r with given ro=5 pm. In spite of a similar shape, there is a substantial 
difference which is due to approximation (7) and neglecting the second order terms 
o[[for) 2 ] in (10]. (In general, p is a linear function of the ratio r/ro.] Left panel quantifies 
the difference between two slopes as a function of r. 


3. Examples for ice crystals 

Is Eq. (3] valid for non-spherical ice crystals? To check this we used the single scattering 
albedo for the Moderate Resolution Imaging Spectroradiometer (MODIS] ice phase function 
available at http:// www . sse c .wi s c . edu / ~ baum/Cirrus /M ODI S V2mode ls.ht m l (see Baum 
et al., 2005a and b] as well as for different types of ice crystals: bullet rosette, plates and 
hollow column provided by Dr. Ping Yang (see Yang et al., 2000]. The results for the MODIS 
ice phase function are illustrated on Fig. 3. Though the correlation coefficients R are not as 
good as for the case of spheres illustrated in Fig. 1, the single scattering albedo for the ice 
crystals clearly shows the same spectral invariant behavior as the one described by Eq. (3]. 
The slope p has a shape similar to Eq. (11], and can be well approximated by 
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where D and Do are effective diameters. (We are not able yet to explain the factor 1.5 
there.] 


4. More Wavelengths 

As mentioned earlier, Eq. (3] is valid for weakly absorbing wavelengths only. The weakly 
absorbing wavelengths are defined here as ones with bulk absorption coefficient fo < 100 
cm 1 (thus for r=10 pm, for < 0.1 and [for) 2 < 0.01]. Left panel in Fig. 4 illustrates the ratio 
between two single scattering albedo spectra o)oA(r]/o)oA(ro] for all wavelengths from 0.2 to 
4 pm and droplet effective radii: r= 20 pm, ro=10 pm. In this plot wavelengths have been 
divided into two spectral regions: 0.2 to 2.5 pm (black dots] and 2.5 to 4 pm (grey dots]. 

For the first region fo< 100 cm' 1 while for the second one fo> 100 cm' 1 (see the right panel 
in Fig. 4 adapted from Fig. 2.25 of Bohren and Clothiaux, 2006]. The ratio for the weakly 



absorbing spectral interval (0.2 - 2.5 pm) was fit with a linear function because 
contributions from the second order terms o((k p r) 2 ) in (7) and (10) can be neglected there. 
The slope of the linear fit p=0.47 is close to the one (p=0.5) determined by Eq. (11) with 
r=20 pm and ro=10 pm. We have also added the four wavelengths used in Fig. 1 (red dots). 
Note that A=3.75 pm is just on the edge between the two spectral regions since /C 3 . 75 s 100 
cm 4 . 


5. Concluding Remarks 

It was shown numerically and confirmed theoretically that the spectra of single scattering 
albedo coox(r) of water droplets of effective radus r at weakly absorbing wavelengths A can 
be well approximated by another spectra cooA(ro) of effective radius ro using Eq. (15) 
assuming that the slope p(r 7 ro) between the ratio cuoA(r)/o)oA(ro) and cooA(r) in the linear 
relationship (3) is known. The slope p(r,ro) can be approximated by a linear function of 
r/ro (see Eq. (11)). In addition to spherical water droplets, it was also shown using the 
MODIS ice phase function (Baum et al., 2005a) that the linear relationship (3) (and thus Eq. 
(15)) is valid for irregular ice crystals of different shapes (Yang et al., 2000). By weakly 
absorbing we mean all wavelengths with bulk absorption coefficient fa < 100 cm' 1 (Bohren 
and Clothiaux, 2006, pg. 113). 

We note that a spectral invariant relationship for single scattering albedo of green leaves 
was theoretically predicted by (Lewis and Disney, 2007). It was found that in the spectral 
interval between 710 nm and 790 nm, spectra of green leaves from different species 
(hazelnut, aspen, jack pine, etc.) were related to one fixed spectrum, called a reference 
spectrum, via Eq. (13); this property was used in classification of forest types from 
hyperspectral data (Schull et al., 2011). 

Finally, linear relationship (3) and its counterpart (15) allow us to relate radiances as a 
function of single scattering albedo cooaM to radiances at other single scattering albedos 
o)oA(r) i n the spectrally invariant relationship between radiances discussed in Marshak et 
al. (2011). This study is yet in progress. 
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Figures 



Figure 1. The ratio of a)oA(r)/oi>oA(ro) plotted against o)oa(c) for four wavelengths A=0. 86, 1.65, 2.13 
and 3.75. Single scattering albedos a)ox(r) are calculated using Mie theory. Droplet effective radius 
r=10, 20, 30, and 40 pm; ro=5 pm. Droplet sizes are assumed to follow a gamma distribution with 
effective variance v =0.1 pm. 



Figure 2. Two slopes p and p app r as a function of r (left panel) and as a scatter plot (right panel). | 

Slope p is from the Mie calculations using four wavelengths (0.86, 1.65, 2.13 and 3.75 pm) while 1 

slope p ap pr is from Eq. (11) where ro=5 pm. | 
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Figure 3. (left) The same as in Fig. 1 but for the MODIS ice phase function (Baum et al., 2005a) with 
effective diameter Do=10 pm. (right) The same as in Fig. 2 (left panel) but for the MODIS ice phase 
function. The approximated slope is described by Eq. (16). 
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Figure 4. (left) The same as in Fig. 1 but for wavelengths A between 0.2 pm and 4 pm with 10 nm 
spectral resolution; r- 20 pm, ro=10 pm. Grey open dots correspond to wavelengths between 2.5 
and 4 pm while black filled dots to wavelengths between 0.2 and 2.5 pm. The linear fit is for the 
black dots only. Red dots correspond to the four wavelengths (0.86, 1.65, 2.13 and 3.75) used in 
Fig. 1. (right) The single scattering albedo o)ox(r) for r=20 pm and the bulk absorption coefficient as 
a function of wavelength. The bulk absorption coefficient, kx, was adapted from Fig. 2.25 of Bohren 
and Clothiaux (2006). The values of kx for liquid water were computed from the refractive index 
and wavelengths compiled by Querry et al. (1991). The red dots correspond to the same 
wavelengths as in left panel. The vertical dash line at 2.5 pm separates two spectral regions. 





